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Abstract. The young stellar object MWC297 is an embedded B1.5Ve star exhibiting strong hydrogen emission lines and a 
strong near-infrared continuum excess. This object has been observed with the VLT interferometer equipped with the AMBER 
instrument during its first commissioning run. VLTI/AMBER is currently the only near infrared interferometer which can 
observe spectrally dispersed visibilities. MWC 297 has been spatially resolved in the continuum with a visibility of 0.50^[j-^Q as 
well as in the Bry emission line where the visibility decrease to a lower value of 0.33 ± 0.06. This change in the visibility with 
the wavelength can be interpreted by the presence of an optically thick disk responsible for the visibility in the continuum and of 
a stellar wind traced by the Bry emission line and whose apparent size is 40% larger We validate this interpretation by building 
a model of the stellar environment that combines a geometrically thin, optically thick accretion disk model consisting of gas and 
dust, and a latitude-dependent stellar wind outflowing above the disk surface. The continuum emission and visibilities obtained 
from this model are fully consistent with the interferometric AMBER data. They agree also with existing optical, near-infrared 
spectra and other broad-band near-infrared interferometric visibilities. We also reproduce the shape of the visibilities in the 
Bry line as well as the profile of this line obtained at an higher spectral resolution with the VLT/ISAAC spectrograph, and 
those of the Ha and Wp lines. The disk and wind models yield a consistent inclination of the system of approximately 20°. A 
picture emerges in which MWC 297 is sunounded by an equatorial flat disk that is possibly still accreting and an outflowing 
wind which has a much higher velocity in the polar region than at the equator The VLTI/AMBER unique capability to measure 
spectral visibilities therefore allows us for the first time to compare the apparent geometry of a wind with the disk structure in 
a young stellar system. 



Key words. Stars: pre-main-sequence, early-type, emission-line: Be stars, individual: MWC 297, planetary systems: proto- 
planetary disks - Infrared: stars - Accretion, accretion disks - technique: interferometric 



1. Introduction 



Table 1. AMBER observation log from 31 May 2004 (UT). 



Pre-main sequence stars in the intermediate mass range, called 
Herbig Ae and Be stars (HAeBe), are observed to be sur- 
rounded by circumstellar material which reveals itself by dis- 
crete emission lines and by continuous excess emission in the 
spectral energy distribution (SED). The spatial distribution of 
this material however has been subject to debate: both geo- 
metrically flat disk models and spherically symmetric enve- 
lope models can reproduce the observed SED. Using viscous 
accretion disk models, Hillenbrand et al. (1992) proposed a 
disk-like geometry for the circumstellar material that however 
needs a central cavity to fit the near-infrared (NIR) excess. 
Considerations on the physical reality of these central cavities 
(Kenyon et al. 1993) and successful fits to HAeBe SEDs by 
dust envelopes (Berrilli et al. 1992; di Francesco et al. 1994) or 
composite envelope-disk models (Miroshnichenko et al. 1999) 
introduced the controversy on the geometry of the circumstellar 
material. 

A break-through occurred with high resolution interfero- 
metric observations of HAeBe stars first presented by Millan- 
Gabet et al. (2001) and Tuthill et al. (2001). These observa- 
tions revealed that the geometries observed in the near-infrared 
bands were closer to either ring-like or spherically symmet- 
ric rather than disk-like. New models introducing passive, star- 
irradiated circumstellar disks with puffed up inner rims were 
found to be consistent with both the interferometric and photo- 
metric observations (Natta et al. 2001; Dullemond et al. 2001) 
although no simultaneous fits were performed. In this scenario 
the inner part of the circumstellar disk edge puffs up due to 
direct irradiation by the central star, roughly equivalent to the 
ring-like structures found with interferometry. The location of 
this inner rim is near the dust sublimation radius, and the sub- 
sequent (dust-free) inner holes are much larger than the ad-hoc 
inner holes proposed in Hillenbrand et al. (1992). Recent stud- 
ies show that there seems to exist a correlation between the 
distance of the inner rim and the luminosity of the central star 
(Monnier & Millan-Gabet 2002). This relation holds well for 
the Herbig Ae and late type Be stars, but breaks down for the 
most luminous HBe stars (Eisner et al. 2004; Monnier et al. 
2005). The latter authors in fact demonstrate that the early-type 
HBe have inner rims that are too close to the star with respect to 
the star's luminosity. The HBe seem to better match the "clas- 
sical" viscous accretion disk model similar to the one initially 
applied by Hillenbrand et al. (1992). A marked diff'erence how- 
ever is that now the inner hole should be filled with optically 
thick gas to effectively shield the dust from destruction at rela- 
tively small distances from the inner rim, assuming that dust is 
what is observed in interferometry. Additional strong evidence 
for a different disk character near the early-type HBe stars is 
found in the spectro-polarimetric fingerprint of Ha lines (Vink 
et al. 2002). These authors prefer to draw the analogy with the 
geometrically flat gaseous disks present near the classical Be 
stars. 
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The geometry of circumstellar material near HAeBe stars 
thus seems to differ between the early-type and late-type mem- 
bers of the group, which is not surprising given the increasing 
interaction between star and disk for the early type stars. For 
the HAe stars a successful working model exists, while on the 
other hand, a disk structure near the HBe stars and their intri- 
cate star-disk interactions still escape a good understanding. In 
this study we present high spatial resolution, intermediate spec- 
tral resolution interferometric observations of the early-type 
Herbig Be star MWC 297using VLTI/ AMBER. This enigmatic 
star appeared in the original HAeBe list compiled by Herbig 
(1960). The star displays a strong emission hne spectrum and 
the character of the underlying photosphere was revealed in the 
detailed study by Drew et al. (1997) to be B 1.5Ve. Its rather 
well determined stellar parameters and its high NIR luminos- 
ity render this star the perfect target to investigate in detail the 
geometry of the circumstellar material near the early type HBe 
stars. 

The paper is organized as follows. Section 2 presents the 
new observations made with AMBER and ISAAC. In Sect. 3, 
we describe the AMBER data processing to produce a reliable 
result presented in Sect. 4 for the visibilities both in the con- 
tinuum and in the Bry line. In Sect. 5 we present an attempt to 
model the environment of MWC 297 with a disk and a wind. 
This model and the consequences are discussed in Sect. 6 and 
summarized in Sect. 7. 

2. Observations 

2.1. AMBER observations 

MWC 297 was observed on 3 1 May 2004 during the second 
night of the first commissioning run of the AMBER instrument 
on the UT2-UT3 (47 m) baseline of the Very Large Telescope 
Interferometer (VLTI). AMBER is the VLTI beam combiner 
operating in the near-infrared (Petrov et al. 2003). The instru- 
ment is based on spatial filtering with fibers and multiaxial 
fringe coding (i.e. AMBER is combining the beams at an an- 
gle which results in fringes modulated in the spatial direction). 
The interferometric beam is anamorphized perpendicularly to 
the fringe coding in order to be injected into the slit of a spec- 
trograph. The instrument can operate at spectral resolutions up 
to 10,000 and efficiently deliver spectrally dispersed visibili- 
ties. 

MWC 297 was measured in the [1980,2230nm] spectral 
range in the MR-K spectral mode (spectral resolution of 1500) 
with 2 elementary detector integration times (DIT) of 3 1 ms 
and 107 ms. Table 1 gives the log of the AMBER observa- 
tions. HD 177756, whose spectral type is B 9V, was used to cal- 
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ibrate the visibilities. Its diameter is 0.60 + 0.06 as computed by 
the ASPRO searchCalib' tool developed at the Jean-Marie 
Mariotti Center. 

The observations of MWC 297 were carried out under spe- 
cific conditions since it was the first AMBER commissioning 
run on the VLTI. The bright Bry line of MWC 297 has been 
originally observed in order to perform a spectral calibration 
of AMBER. Detailed analysis of the commissioning data from 
this run and later ones has shown that the optical train of the 
UT telescopes were affected by non-stationary high-amplitude 
vibrations. Because of the small number of observations on 
MWC 297, and because the amplitude of the vibrations might 
undergo rapid variations, the calibration of our measurements 
must therefore be regarded with care. In order to investigate 
the errors on the visibility measurements and check their con- 
sistency, we have used different data reduction methods and 
different data selection schemes which are described below. 
However, we would like to point out that the vibrations do not 
impact the spectral dependence of the visibility since they af- 
fect the spectral range as a whole. 

2.2. ISAAC observations 

MWC 297 was observed in service mode on 13 July 2004 
with the ESO VLT UTl telescope under modest seeing con- 
ditions (~ 1.5" in the visual). The ISAAC near-infrared spec- 
trograph was employed in the short wavelength medium resolu- 
tion mode with a 0.3" wide slit. This instrument setup delivered 
a resolution of A/AA ~ 8900 at the Bry wavelength. The raw 
data were flat-fielded, wavelength-calibrated and corrected for 
telluric absorption using standard techniques and observations 
taken from the ISAAC calibration plan by the ESO staff as- 
tronomers. A detailed account of the observations and the data 
reduction is given in Garcia Lopez et al. (2005). 

2.3. Existing photometric and interferometric data sets 

Extensive photometric and interferometric data exist for 
MWC 297. Broad-band photometric data were collected and 
presented by Pezzuto et al. (1997). This data set consists of 
UBVRI from Bergner et al. (1988), JHKLMN from Berrifli 
et al. (1992), and Q-band data from Simon (1974). In the 
mm/submm wavelength regime (0.35-1. 3mm) the dust contin- 
uum measurements are taken from Mannings (1994). At ra- 
dio wavelengths (6 cm), MWC 297 has been observed by Drew 
et al. (1997). Existing NIR interferometric data for MWC 297 
consist of two sets. IOTA H-band continuum data were pre- 
sented by Miflan-Gabet et al. (2001), and PTI K-band contin- 
uum data were published by Eisner et al. (2004, upper limits 
only). 

3. AMBER data processing 

3.1. Raw data 

AMBER follows the standard data flow system implemented 
at ESO/VLT. During data acquisition, the software records the 

' ASPRO is available at http : //mariotti . £r 




Fig. 1. Short-exposure image of the MWC 297 signal on the AMBER 
detector. X-axis corresponds to the spatial extension of the beams and 
Y-axis to the wavelength. First column (Dk) corresponds to the dark, 
the second (PI) and third (P2) ones are the beams from the first and 
second telescope resp., and finally the last column (IF) shows the 
fringes obtained by superposition of the two beams. The bright row 
is the Bry line at 2165.6 nm. 

images of the spectrally dispersed fringes as well as those of 
the telescope beams. Figure 1 shows a typical image retrieved 
from the detector subsystem-. From the left to the right, the 
first column contains an estimation of the dark, the second and 
third ones are the spectrophotometric measurements from each 
telescope in order to evaluate their respective contribution to 
the interferogram, and the last one is the spectrally dispersed 
fringe pattern, corresponding to the interference between the 
two telescope beams. Each exposure (or file) consists of one 
hundred of these elementary frames. 

One can notice the bright line near the top of the figure 
and the dark lines at the bottom. They correspond to the Bry 
emission line at 2. 1656 /^m and to the telluric absorption lines 
around 2.05 and 2.08 fim respectively. 

3.2. Data reduction 

Since we are still in an early stage of AMBER observations, 
we decided to use two different methods to extract the raw vis- 
ibilites. We present and discuss them briefly. 

- The standard data reduction method developed and opti- 
mized for AMBER is called P2 VM for Pixel-To- VisibiHties 
Matrix (Millour et al. 2004; TatulH et al. 2005). The 
P2VM is a linear matrix method which derives raw visi- 
bilities from AMBER data for each spectral channel. The 
P2VM is computed after an internal calibration procedure 
which is performed every time the instrument configura- 

- As a matter of fact, this image is not exactly how it appears on the 
detector, since the four columns come from more widely separated 
regions of the detector. 
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tion changes. The complex coherent fluxes^ are given by 
the product of the fluxes measured on each pixel of the de- 
tector by this P2VM. We estimate the visibility for each 
spectral channel by integrating the square amplitude of the 
complex coherent fluxes over a certain number of frames 
and then by normalizing it by the total flux obtained for the 
same set of frames. The number of frames used can vary 
between 1 and the total number of frames in one exposure. 
This parameter is called hereafter the binning. 

- A Fourier transform (FT) technique (see Ohnaka et al. 
2003, appendix A) has been used. The advantage is that 
it does not depend on the internal calibration of the instru- 
ment, described above. To enhance the signal, we applied 
a sliding average of seven spectral channels which has no 
effect on the continuum visibilities, but allows us to retrieve 
only upper limits of the visibility in the Bry line. The visi- 
bility is computed using power spectral densities averaged 
over a variable number of frames. 

Both methods compute squared visibilities and therefore an im- 
portant step is to subtract the visibility bias, i.e. the part of the 
visibility due to photon and read-out noises. Because of the 
presence of vibrations, we are not completely sure that this bias 
is perfectly estimated and therefore we use very conservative 
errors. 

3.3. Data selection 

Because of the specific conditions of our observations, one crit- 
ical step in the processing is the selection of the best data within 
the whole set. We based our selection of individual frames or 
exposures on the value of the visibility signal-to-noise ratio 
(SNR). 

We used two different methods of selection based either on 
exposure or frame selection. 

- The selection over the exposures consists of selecting a cer- 
tain fraction of the best exposures. The visibilities for each 
exposure are computed using all frames. The final value is 
the average of the visibilities over the selected exposures 
weighted by the statistical errors. 

- The selection over the frames consists of selecting a certain 
fraction of the best frames within each exposures acquired 
on the object. The final value of the visibility is the average 
of the visibilities of all selected frames weighted by their 
errors. 

The exposure selection has been used only with the P2VM 
method whereas the frame selection has been used with both 
reduction methods. In all cases, the selection has been applied 
to the object and its calibrator and the resulting calibrated visi- 
bilities are summarized in Table 2. 

The threshold for the selection process is an important pa- 
rameter. Selecting with a strong criteria (e.g. < 5%) would not 
provide reliable statistics for the biases on the square visibil- 
ity. Selection with a soft criteria (e.g above > 50%) may re- 
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Fig. 2. Spectral dependence of the visibility as measured with 
AMBER for MWC 297 around the Bry line. 

tain influences by a reduced fringe contrast due to telescope vi- 
brations. Therefore we chose 10% for the frame selection, and 
25% (resp. 35%) for the exposure selection of the 31 ms (resp. 
107 ms) data set because of the small number of exposures (see 
log of observations). 

4. Results 

4.1. Continuum visibilities 

The continuum visibiUties computed with the different algo- 
rithms and selection schemes are summarized in the right col- 
umn of Table 2. We find discrepancies between the two meth- 
ods that we are not able to explain in this very early stage 
of the VLTI and AMBER. These discrepancies might be re- 
lated to the presence of vibrations in the VLT UT coude trains 
as mentioned in Sect. 2, or to remaining imperfection in the 
estimation or subtraction of the biases. We propose to take 
the following continuum visibility with a relatively large er- 



ror: Vcon, = 0.50 



+0.08 



This confirms that the environment of 



Coherent flux is the degree of complex coherence not yet normal- 
ized by the flux. 



^cont - u.-yu_Q jo. 

MWC 297 is spatially resolved at the level of a few milliarc- 
seconds (Eisner et al. 2004). The uniform disk diameter cor- 
responding to the visibility measured is 7.0 + 0.9 mas, corre- 
sponding to 1.75 + 0.23 AU at 250 pc. 

4.2. Spectral variation of the visibilities 

All methods and selection schemes give consistent visibilities 
for the Bry emission line at 2.1656yL(m: V\ine - 0.33 + 0.06 
(see Table 2). Fig. 2 shows the variation of the visibility with 
wavelength for the 3 1 ms data. For this figure, the continuum 
and line visibilities were calculated using the P2VM and the 
best exposure selection criterion. 
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Fig. 3. Comparison of Bry observed with AMBER (solid line) and 
ISAAC (dotted line). The dashed line corresponds to the ISAAC spec- 
trum convolved by the response function of 
AMBER. 



Fig. 4. Sketch of the model including an optically thick disk and an 
outflowing wind (edge-on view). The receding part of the wind is only 
partly visible because of the screen made by the optically thick disk. 



The visibility from the pure Bry emission must be corrected 
from the influence of the continuum by the following relation- 
ship: 

^line^line ^cont Vcont 

Vbit = — (1) 

line cont 

where Fbij = f'line - ^^cont- Since the ratio of the line flux to the 
continuum flux is Fiine/^^cont - 2.2, the visibility of the region 
emitting the Bry line is Vbtj = 0.19 + 0.03, corresponding to 
a uniform disk diameter of 9.9 + 0.3 mas or 2.5 + 0.08 AU at 
250 pc. Therefore, the size of the Bry emitting region is 40% 
larger (in terms of uniform disk diameter) than the size of the 
region contributing to the continuum. 

4.3. ISAAC spectrum 

The AMBER spectrum at medium resolution CR » 1500) con- 
tains the Bry line. Its profile has been compared to the ISAAC 
high spectral resolution CR x 8900) spectra obtained at nearly 
the same epoch in May and luly 2004 respectively. A compar- 
ison between the two lines is given in Fig. 3. The equivalent 
width of the AMBER spectrum is 2.14nm and is compatible 
with the 1 .58 nm measured from the ISAAC spectrum given the 



spectral variability of MWC 297. We note that the resolved in- 
frared Bry emission lines presented in Murdoch & Drew (1994) 
are symmetric rather than double-peaked. 

The ISAAC spectral resolution is larger than the one 
achieved with AMBER. Nevertheless, in order to compare our 
two data sets the ISAAC line profile was convolved with a 
Gaussian function corresponding to the AMBER spectral reso- 
lution. The shape of both profiles are quite similar (see Fig. 3) 
but since the higher spectral resolution of the ISAAC spectrum 
provides more details on the kinematics within the circumstel- 
lar envelope, we have used, in the following, the ISAAC Bry 
line profile in order to constrain the outflowing wind model 
(Sect. 5.2). 

5. Modeling 

In this section we present the modeling of the large body of 
interferometric, spectroscopic and photometric data that exists 
for MWC 297. The modeling is done by applying two differ- 
ent codes, one for an optically thick disk and one for a stellar 
wind. The disk code is designed to model the continuum ra- 
diation, whereas the stellar wind code reproduces the strong 
emission lines. Figure 4 represents a sketch of the combined 
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Table 3. Thie best-fit accretion disk model parameters found by simul- 
taneous fitting of SED and visibilities. 



Accretion rate (Mace) 
Inner radius (^i„) 
Outer radius (^oui) 
Inclination (;) 
Position angle (0) 



- 1 X 10"^ Mo yr-' 
0.5 + O.lAU 
55±5AU 
15 ±5° 
56 ±7° 



model, where the optically thick disk and the outflowing wind 
are spatially independent. 

5.1. Continuum radiation: optically thick disk 

The disk model is based on the simple model already used by 
Malbet & Bertout (1995) and Malbet et al. (2005). This disk 
model consists of an axisymmetric radial analytic disk struc- 
ture which is heated both by stationary accretion and stellar 
irradiation. The disk is in hydrostatic equilibrium and non self- 
gravitating. The accretion flux is following the standard power 
law for a viscous disk (Lynden-Bell & Pringle 1974; Shakura & 
Sunyaev 1973). The emitted continuum flux is produced by the 
emission of optically thick but geometrically thin black-body 
radiating rings. It produces an SED, and its spatial distribution 
can be Fourier transformed, which reduces to a Hankel trans- 
form for the radial symmetry of a disk, to obtain interferometric 
visibilities. 

The model effectively has five free parameters: the accre- 
tion rate Mace, the inclination ;, the position angle 0, and the 
distances of the inner and outer edge of the disk, Rin and Rom re- 
spectively. Modeling is strongly constrained by the accurately 
determined properties of the central star by Drew et al. (1997): 
a stellar radius of 6.12 Rq, a mass of 10 M©, an effective temper- 
ature of 23, 700 K for a distance of 250 pc with a visual extinc- 
tion of Ay = 8 mag. The outer edge of the disk is strongly con- 
strained by the mm/submm data point: Mannings (1994) finds 
an outer disk radius of ^ 60 AU. In short, simultaneous fitting 
of the SED and the visibilities lead to balancing out the mass 
accretion rate, inner disk radius and the inclination. In practice, 
this is done iteratively by first fitting the SED, with / and 9 set 
to zero. We deem the R- and /-band continuum measurements 
not reliable given the very strong Ha emission; these two are 
therefore not taken into account in the fit. In the first step we 
obtain Mace and R^^- These are then used in a separate fit of the 
visibiUties as function of / and 6. The results are used again for 
the first step and this procedure is repeated until convergence. 
In Table 3 we list the obtained best-fit model values. The mod- 
eled SED and visibilities are found to be in good agreement 
with the observed values, as shown in Fig 5. 

We probed the sensitivity of these fits by varying the central 
star parameters, according to the uncertainties given by Drew 
et al. (1997). They derived half a spectral subtype uncertainty, 
and a distance error of 50 pc. Fitting the SED using the extreme 
values for the central star, we find the uncertainties quoted in 
Table 3. Especially the mass accretion rate is far from well de- 



termined. If the central star would be of type B2 at a distance of 
200 pc, the required mass accretion rate is only ~ 10"^ M© yr"'. 

5.2. Emission lines: optically thin outflowing wind 

In our model, the emission lines are produced in a circumstellar 
gas envelope. In order to model this line profile and the corre- 
sponding visibilities, we have used the SIMECA code (Stee & 
de Araujo 1994; Stee et al. 1995). This code computes clas- 
sical observables, i.e. spectroscopic and photometric ones but 
also intensity maps in Balmer lines and in the continuum in 
order to obtain theoretical visibility curves. The main assump- 
tions are that (i) the envelope is axisymmetric with respect to 
the rotational axis, (ii) no meridian circulation is allowed, (iii) 
the physics of the polar regions is well represented by a CAK 
type stellar wind model (Castor et al. 1975). The solutions for 
all stellar latitudes are obtained by introducing a parametrized 
model (power of sinus function) constrained by the spectrally 
resolved interferometric data. Depending on the value of the 
chosen terminal velocity at the equator, the equatorial region 
can be dominated either by Keplerian rotation or by expansion. 

Since the SIMECA code has originally been developed to 
model the circumstellar environment of classical Be stars, we 
had to modify the code in order to interface SIMECA with 
the optically thick disk model described previously. We have 
implemented three changes (the equations describing the wind 
model are recalled in Appendix A): 

1. The wind is no longer computed from the equator to the 
pole, but the computation occurs in a bipolar cone defined 
by a minimal angle allowing the disk to be present (see 
sketch in Fig. 4). The disk model tells us that the opening 
angle is between 1.8° at the inner radius and up to 3.2° at 
the outer radius. We used a minimum angle of 4°. Therefore 
the equatorial terminal velocity corresponds to the termi- 
nal velocity at this minimal angle from the equatorial plane 
at the interface between the accretion disk and the stellar 
wind. 

2. The disk hides the receding part of the wind. In Fig. 4, the 
part of the wind which is not visible from the observer is 
not taken into account in the outgoing flux. 

3. Although the disk emission contributes less than 1% com- 
pared to the star flux in the visible (i.e. also in the Ha and 
H/3 lines) and can be neglected, at 2. 1656 fim the disk emis- 
sion is 6.4 times larger than the stellar flux. This contribu- 
tion is taken into account in addition to the excess due to 
free-free emission from the outflowing gas when comput- 
ing continuum normalized intensities and also for the com- 
putation of the visibilities. 

To be consistent with the disk model (although it might not 
be fully physical), we also use the same temperature and rota- 
tional velocity laws, respectively T{r) oc r^^^'^ and a Keplerian 
rotation. 

After running hundreds of simulations in order to constrain 
the physical parameters of the wind and to test the geometri- 
cal and kinematic hypothesis described later in Sect. 6.2 , we 
find a successful simultaneous fit to the ISAAC Bry line pro- 
file (see Fig. 6) as well as Ha, H/3 profiles compatible with 




Wavelength (^m) Baseline (m) 

Fig. 5. Result from the optically thick disk model. Left panel: observed and modeled SED for MWC297. The full dots are the continuum 
measurements taken from Pezzuto et al. (1997), also included are the ISO SWS/LWS spectra. Dotted line is the star, dashed line the accretion 
disk, and the full line the resulting total flux of the model. Right panel: resulting best-fit model radial squared visibilities compared with 
continuum squared visibilities observed with AMBER, IOTA and PTI. Full line and IOTA data are in the //-band, dashed line and AMBER/PTI 
are in the A'-band. PTI values are upper limits. 
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Fig. 6. Double peaked Bry profile observed by ISAAC (full line) and 
modeled with the outflowing wind model (dotted line). We have also 

plotted the cumulative contribution of the accretion disk (dashed line) ''■ visibility observed with AMBER (points with error bars) 



and of the star (dashed-dotted). 



and the one obtained from the outflowing wind model (full line). 



Table 4. The best-fit model parameters for the outflowing wind model. 
Most parameters are self-explanatory (see Appendix A for details). C\ 
is the ratio between the polar and equatorial mass flux; m\ is the ex- 
ponent of the mass flux law as function of latitude; m2 is the exponent 
of the latitude dependent terminal velocity law. 



Photospheric density 

Equatorial rotational velocity 

Polar terminal velocity 

Terminal velocity above disk 

Polar mass flux 

Ci 

mi 

nil 

Inclination angle (/) 



1 ±0.5 X 10'- cm- 



400 
600 
70 ± 
3.2 H 



t 50 km s 
t 50 km s 
20kms-' 
:0.2x 10 



'Moyr- 



0.25 ± 0.05 
30 ± 10 
10 ±2 

25 ±5° 



Drew's observations. The best-fit model parameters are given 
in Table 4. The outflowing wind model reproduces successfully 
the AMBER measured drop in visibility across the Bry line as 
shown in Fig. 7. 

The computed Wa and Wfi intensities are respectively 130 
and 11.5 compared to 120 and 11 obtained by Drew et al. 
(1997). This corresponds to an intensity ratio ly^allyip of 11.3 
compatible with the Drew's ratio of 11 regarding the stel- 
lar variability and the non-simultaneity of the measurements. 
Moreover we were able to reproduce quite well the shape of 
these lines (i.e width) and the double peaked emission of the 
Bry line. The asymmetry of the two peaks of the Bry line is 
also successfully reproduced thanks to the introduction in the 
SIMECA code of the opacity of the disk (point 2 of SIMECA 
modifications). Nevertheless the agreement is not perfect in the 
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red wing of the profile probably due to our ad-hoc way of in- 
terfacing of the wind and the disk. 

The global shape of these lines and their intensities are very 
sensitive to the inclination angle. We were able to determine the 
inclination of 25 + 5° from the fitting of the line profiles with 
the SIMECA code. We note that the wind parameters given in 
Table 4 are different from classical Be stars. The mass flux ratio 
(Ci) is 0.25 whereas for Be stars the allowed range is between 
10 and 100. Knowing that the density in the envelope is propor- 
tional to the mass flux and inversely proportional to the radial 
velocity we obtain that the gas density in the polar region is 
twice larger than the one at the interface between the accretion 
disk and the stellar wind. The highest mass loss occurs along 
the polar direction. 

6. Discussion 

We have presented extensive existing data and new AMBER 
and ISAAC data on the Herbig Be star MWC 297. We have si- 
multaneously modeled interferometric and non-interferometric 
data. The continuum emission is generated by an optically thick 
disk heated by stellar irradiation and accretion, whereas optical 
and NIR emission lines find their origin in an outflowing stellar 
wind. 

The modeling presented in the previous section, although 
rather successful, raises new questions on the physics of the 
circumstellar environment of intermediate-mass young stars. In 
this section, we discuss them by addressing first separate issues 
about the disk and the wind, and then those raised by the com- 
bined model. 

6.1. Physical conditions in the disk 

Our continuum observations placed in the context of young 
stellar objects allow us to bring new elements into the discus- 
sion of protoplanetary disk models. 

6.1.1. Reprocessing versus accretion in MWC 297 

The mass accretion rate is hardly constrained by the viscous 
accretion disk model. In the first place because the viscosity in 
the accretion disk is quite poorly understood, and therefore the 
quoted accretion rate cannot be assessed with real certainty. In 
the second place, we obtain two orders of magnitude variation 
in the value of the accretion rate, by varying the stellar param- 
eters according to their uncertainties (see Sect. 5.1). 

The ratio between the reprocessed energy from the stellar 
irradiation and the viscous heating due to the accretion along 
the disk varies from4x down to 3x(Macc/10"^ Mq yr"')"' from 
the inner part to the outer part of the disk. Therefore for accre- 
tion rates lower than lO^^Moyr"', the reprocessing energy is 
always dominant compared to the accretion energy. In fact the 
transition between the two types of flux in the disk lies near a 
Mace - 3 X 10"^ Mq yr"', for the most probable values of the 
stellar parameters. 

This means that the disk is predominantly a reprocessing 
disk. However as pointed out by Lachaume et al. (2003, see 



their Fig. 4), the emerging flux might be dominated by the stel- 
lar heating and the vertical structure by the viscous heating at 
least for the inner radii (less than a few AUs). 

6.1 .2. Inner radius of the ciisl< 

The best fit of our disk model is found with an inner disk 
boundary at 0.5 AU. This result has to be compared to the 
analysis of Millan-Gabet et al. (2001). With an assumption of 
d - 450 pc, these authors found a ring diameter of about 10 AU 
and an inner disk radius of 0.9 AU. Our result is therefore con- 
sistent with their disk model since we have adopted the distance 
d = 250 pc of Drew et al. (1997). 

As already pointed out by other authors (Eisner et al. 
2004; Monnier et al. 2005), the inferred inner disk radius for 
MWC 297 is not compatible with the dust evaporation distance 
from the star. Using the disk model with the "puff'ed up" inner 
rim (Dullemond et al. 2001; Isella & Natta 2005), the inner ra- 
dius of the disk would be located at about 3 AU from the star, 
even in the most favourable hypothesis that the dust evapora- 
tion temperature is of 2000 K. Moreover, MWC 297 seems not 
to be a special case, since this behavior is common to the other 
two early HBe stars (Z CMa A with a spectral type B and 
V 1685 Cyg with spectral type B 3) observed with near infrared 
interferometers (Monnier et al. 2005). One possible explana- 
tion may be that the gas is the dominant source of absorption 
in the inner part of the disk, preventing dust grains from evap- 
oration near to the star 

Using Rosseland gray opacities (Ferguson et al. 2005), we 
checked that the disk is always optically thick both for an ac- 
creting disk (accretion rate higher than lO^^Moyr"') or for a 
reprocessed disk (disk density at the disk edge higher than 1% 
of the typical value of the solar nebula'*; Hayashi 1981; Davis 
2005). 

Even with the last assumption, the structure of the inner 
region of early HBe stars is not totally understood. The physical 
reason why the optically thick gas disk should be truncated at 
0.5 AU from the star is not clear. 

6.1 .3. Ionization of the disk 

We estimate if the circumstellar disk is susceptible to magne- 
tohydrodynamical (MHD) instabilities, and thus if the disk is 
in active accretion state by means of magneto-rotational insta- 
bilities (MRI; Balbus & Hawley 1991). Given the large inner 
radius of 0.5 AU, it is not expected that the disk is accreting, 
otherwise the region interior to the disk inner radius would 
quickly be filled. Indeed, in the previously discussed hydro- 
dynamic case, we saw that the disk seems to be dominated by 
reprocessing. 

Whether MHD instabilities in a given circumstellar disk 
operate or are suppressed can be estimated with the mag- 
netic Reynolds number (see e.g. Gammie 1996). By setting 
the Reynolds number equal to one, a threshold disk ionization 
fraction, of 3 x lO"'-' is found. In this computation we 

adopted the same parameters for the MWC 297 star-disk sys- 

Solar nebula density law used: Es A. = 1.7x10^ (r/lAU)"'-^ g/cm^ 
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tem as in the previous sections. The ionization fraction thresh- 
old is essentially determined by the ionization rate and the 
particle density in the disk. Considering a high (IO^^Mq yr ') 
and a low (10 "^Mgyr"') mass accretion rate for MWC 297, 
we can use Fig. 1 of Inutsuka & Sano (2005) to deduce that 
the required ionization rate in a dusty disk should range be- 
tween 10"'"* and 10"'^ s"' or between 10"'*' and 10"" s"' re- 
spectively. Given that the cosmic ray ionization rate is close to 
10^'^ s"\ we conclude that cosmic rays alone are insufficient to 
incite and/or sustain the required fractional ionization for MRI 
to operate. 

Since ultraviolet radiation is generally considered to be in- 
efficient due to the very small attenuation length, X-ray ioniza- 
tion may play an important role in the ionization structure of 
the disk. The very strong X-ray activity reported for MWC 297 
with ASCA (Hamaguchi et al. 2000) should be interpreted 
with care given the reasonably strong case for source confusion 
(Vink et al. 2005). An inner disk radius of 0.5 AU is thus prob- 
ably not susceptible to MHD instabilities. We conclude thus 
that from a MHD and a hydrodynamic point of view an inert 
reprocessing disk may exist. We note however that if in the past 
the disk of MWC 297 was in fact a magnetically active disk in 
which MRI operated, that then the concomitant MRI turbulence 
could have sustained the ionization degree above the threshold 
level even in the presence of dust (Inutsuka & Sano 2005). 

6.2. Geometry of the wind 

First, we investigate why the standard Be model of a rotating 
wind is unsuccessful to reproduce the MWC 297 data. In this 
model (Stee & de Araujo 1994; Stee et al. 1995), the Bry emis- 
sion would originate from a rotating flattened envelope with no 
or little radial expansion. We can test the consistency of this 
model by comparing the size of the emitting region given by 
the kinematics revealed by the splitting of the Bry emission 
line into two peaks and the radius of this same region inferred 
by the drop in visibility in this Bry line. 

On the one hand, the projected rotational velocity of the or- 
biting material is given by the separation of the two Bry peaks 
and is proportional to the maximum extension of the emitting 
region (Huang 1972; Hirata & Kogure 1984). Using a rotational 
velocity law of the form v(r) = vo (r//?*)"*', where vo is the star 
rotational velocity at the equator, r is the distance to the star, 
and X is the exponent of the rotational velocity law, the separa- 
tion between the two peaks is proportional to 2vo sin / (r//?*)"*. 
Therefore the maximal radius of emission of the Bry line is 
~ (2vo sin //t/)'^-*', and the stellar parameters given by Drew 
et al. (1997) and a Keplerian rotation yield ~ 1407?*. 

Since the typical radius of emission of the Bry emission 
is ~ 43/?* (see Sect. 4.2), the size of 140/?* derived from the 
Keplerian rotation of a rotating stellar Be wind is not consistent 
with our interferometric measurement. 

In fact, in order to obtain an emitting zone radius compat- 
ible with our interferometric measurement we need to assume 
a rotational velocity of less than 250 km s"' or a rotation law 
far from the Keplerian one. In total, we ran 50 models in our 
attempt to fit the Bry, Ha, Hfi lines and the visibility in Bry 
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Fig. 9. Continuum //-band squared visibilities obtained with the disk 
model compared to the IOTA visibilities observed during the transit of 
the object over the sky with the IOTA S15N35 baseline. 

simultaneously using this type of wind model. None of these 
models was able to fit correctly the Bry line profile displayed 
in Fig. 3 nor the intensity ratio of the three hydrogen lines. In 
addition, the Ha and H/3 emission lines require to form at least 
partially in a strong stellar wind to reproduce the 600 km s"' 
line width in contrast with the Bry line displaying a 60kms"' 
double-peaked profile. This is unlikely in the same environ- 
ment. 

The outflowing stellar wind coupled with an optically thick 
disk gives a better explanation of the kinematical and geomet- 
rical morphology of MWC 297. In Fig. 8, we present the inten- 
sity maps of MWC 297 computed with our model seen edge-on 
in the three hydrogen lines in order to better localize the region 
of emission of the lines. The emission from Ha and H/3 orig- 
inates from a large and somewhat spherical region where the 
velocity can reach up to 600 km s"'. The Bry line is confined 
into a narrower region just above the optically thick disk where 
the velocity is dominated by the disk Keplerian rotation and a 
terminal velocity of 70kms"' (see parameters in Table 4). In 
fact Murdoch & Drew (1994) have already pointed out that the 
Bra/Bry line flux ratio shows an increase at low velocity that 
cannot be interpreted as an outwardly-accelerated wind model. 

6.3. Disk and wind interaction in MWC 297 

Even if the disk and the wind have been modeled separately, 
these two physical phenomena need to be compared at least at 
the parameter level. We find a good correspondence between 
the inclinations deduced from the disk model (/ = 15+5° 
from the visibility fit) and from the wind model (; - 25+5° 
from the asymmetry of the Bry line peaks). The disk inclina- 
tion is strongly constrained by the IOTA visibilities spread out 
over a large range of hour angle as shown in Fig. 9. The wind 
mass loss rate is smaller than the maximum disk accretion rate 
by several orders of magnitude, and therefore compatible with 
most disk/wind theories but also with most observations of out- 
flows and accretion activities. As shown in Sect. 6.1.3, the inner 
radius is hardly ionized and ffius will not contribute to the hy- 
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Fig. 8. Edge-on intensity maps of the wind in the computed Hey, H/3, Bry lines. 
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Fig. 10. Pole-on intensity maps of the wind Bry emission (left panel) and of the A^-band disk continuum emission (center panel). Right panel 
shows a radial cut of these intensity maps with the Bry wind in solid line and the continuum disk in dashed line. 



drogen lines, confirming our choice to model Bry only in the 
wind. 

Figure 10 shows the pole-on intensity maps of the disk 
model in the continuum and of the wind in the Bry line, as well 
as their respective intensity profile. This is a graphical explana- 
tion of the visibilities observed by AMBER: the wind angular 
extension in the Bry line is larger than the disk apparent size, 
and therefore the visibility is smaller within the line. 

Can the result obtained with these observations constrain 
the nature of the wind? We recall that in protoplanetary disks, 
two main classes of disk wind models have been proposed de- 
pending on the geometry of the magnetic field lines. 

- The disk wind initially proposed by Blandford & Payne 
(1982) and studied in detail for example in Casse & Ferreira 
(2000). In this theory the magnetic field is attached to the 
disk and the particles which accrete to the star may be 
launched into a wind eventually collimated at a certain dis- 
tance of the star. 

- The X-wind of Shu et al. (1994) has a different geometry 
and originates at the corotation radius of the disk where all 
magnetic field lines are localized in the equatorial plane. 

With the present state of data, we cannot even distinguish be- 
tween a stellar wind and such disk winds since we are unable to 
recover a precise kinematical map of the wind. We note how- 



ever that our model with Keplerian rotation does not respect the 
conservation of angular momentum and can only be sustained 
if additional MHD energy is injected in the wind. More spec- 
tral resolution with AMBER will help to answer the question, 
especially using the 10, 000 spectral resolution mode. 

In conclusion, we can claim that the models of disk and 
wind are compatible and the combination of them is probably 
close to the reality. A complete and self-consistent modeling 
of the environment is out of the scope of this paper but would 
allow to better constrain the physical processes in action both 
in the disk and the wind, like accretion and ejection, at least 
from the observational point of view. 

6.4. Inclination discrepancy 

The modeling presented in this paper does a reasonably good 
job in reproducing nearly all the observational data and pro- 
duces fiducial physical parameters for the circumstellar envi- 
ronment of MWC 297. However, we derive an inclination of 
~ 20° for the system, which is not consistent with a near edge- 
on orientation as proposed by Drew et al. (1997). The latter is 
inferred from the photospheric lines that indicate a 350 km s"' 
projected rotational velocity. An inclination of 20° would lead 
to a rotation above the break-up velocity. 
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The controversy becomes even more pronounced given the 
Ha spectro-polarimetric results by Oudmaijer & Drew (1999). 
These authors find no effect across the line, a result consistent 
with a symmetric emission zone for Ha. On the contrary, the 
morphology of the radio emission (Drew et al. 1997) favors a 
bipolar interpretation and again would point to a high inclina- 
tion of the system, although not completely inconsistent with 
a symmetric geometry on smaller scales (Oudmaijer & Drew 
1999). 

A possible reconciliation of the seemingly contradictory 
observables and derived physical parameters is a tidally in- 
duced warped geometry for the circumstellar disk (Terquem & 
Bertout 1993) due to a stellar companion located at a distance 
of a few hundred AU. Recently, Vink et al. (2005) reported the 
presence of a close-by late-type companion at a projected dis- 
tance of 850 AU. By converting the //-band magnitude of the 
Vink et al. (2005) companion star to a spectral type A2 V and 
by calculating the steady-state deformation of the MWC 297 
circumstellar disk due to the star at 850 AU, (Eq. 25 of Terquem 
& Bertout 1993), we find a deformation of less than 1/1000 
even at the edge of the disk at 60 AU. Clearly this star would 
not be the right candidate, but we feel that the possibility of a 
tidally deformed disk may not yet be ruled out given the many 
companion stars found near MWC 297 (Testi et al. 1999; Vink 
et al. 2005). An even closer companion can not be excluded, a 
hint of which may be the CHANDRA X-rays found coincident 
with the position of MWC 297. 

7. Summary and perspectives 

We have presented first spatially and spectrally resolved obser- 
vations of the disk/wind interaction in the young stellar system 
MWC 297 with the VLT interferometer equipped with the in- 
strument AMBER. The unique spectral capability of AMBER 
has allowed us to measure for the first time the visibility in 
the Bry emission line in a young stellar object and to compare 
it to the continuum visibility. The continuum visibility in the 
/T-band drops from 0.50 to 0.33 in the Bry emission line of 
MWC 297. The spectrum obtained with AMBER is consistent 
with a double peaked spectrum observed with ISAAC on the 
VLT, where the peaks are roughly separated by approximately 
60kms"'. 

We have successfully modeled the circumstellar environ- 
ment of MWC 297 using an optically thick geometrically thin 
disk and an outflowing stellar radial wind having an increas- 
ing outflowing velocity starting from the surface of the disk up 
to the pole. This combined model is able to reproduce many 
observational features like the continuum visibilities measured 
by AMBER as well as the visibilities in the Bry line, together 
with the shape of the SED over more than three orders of mag- 
nitude of the wavelengths, the broad-band visibilities obtained 
by other infrared interferometers as well as the Ha, li/3, and 
Bry line profiles. 

We have discussed our result in the light of more sophis- 
ticated models. We showed that the inner radius is not deter- 
mined by the dust sublimation distance, and is unlikely ion- 
ized by cosmic rays only. The disk flux is mainly driven by the 
stellar reprocessing although we cannot rule out that the accre- 



tion process may play a role on the vertical structure. We have 
shown that our AMBER observations interpreted by our model 
predict emission of the Ha and H/3 lines in the polar regions 
whereas the Bry emission arises mainly from the region just 
above the surface of the disk. AMBER continuum and line ob- 
servations point both toward an system inclination of approxi- 
matively 20°. 

We are not yet able to constrain the exact nature of the 
wind and the type of connection with the disk, but we expect 
that future observations with AMBER will bring new pieces 
in the understanding of MWC 297. As a matter of fact, the 
vibrations in the VLTI UT coude train should be diagnosed 
and fixed in a short term, the VLTI equipped with a fringe 
tracker, allowing AMBER to be operated at its highest spec- 
tral resolution (~ 10,000) which will give new kinematical 
information on this interesting and intriguing region around 
MWC 297. Another direction of investigation would be to carry 
out AMBER observations with three telescopes, in order to 
measure the closure phase with different baseline configura- 
tions, and therefore measure possible departure from centro 
symmetry of the material around MWC 297. 
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Appendix A: Wind model equations 

Following Stee & de Araujo (1994) and Stee et al. (1995), the 
mass flux is parametrized as: 

O(0) = (Dp„i, + [((D,q - Op^i,) sin"" (6)], (A. 1) 

where mi is the first free parameter which describes the vari- 
ation of the mass flux from the pole to the equator. The ratio 

' ftp : //ftp-icf . llnl . gov/pub/Yorick 
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between the equatorial and polar mass flux is C\ - Oeq/<l>,poie- 
The values of Ci are typically between 10' and lO'* (Lamers & 
Waters 1987). Equation (A.l) can be rewritten as: 



3>(0) = 1'poie[l+(C, -l)sin""(0)]. 
The expansion velocity field is given by: 

v,(r,0) = vM + [vuo) - Voiom - -y, 

r 

with y - 0.86 and 

cD(0) (Dpoie[l+(Ci -l)sin""(^^)] 



Void) 



Po 



(A.2) 



(A.3) 



(A.4) 



The second free parameter m2 is introduced in the expres- 
sion of the terminal velocity as a function of the stellar latitude: 

Voo(0) = Voo(pole) + [VocCeq) - VooCpole)] sin'"^(0). (A.5) 

Finally the density distribution in the envelope is given by the 
equation of mass conservation: 



p{r,0) 



CD(0) 



(A.6) 
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